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We have used a dc Superconducting QUantum Interference Device to measure the magnetic 
flux noise generated by the equilibrium vortex density fluctuations associated with the Kosterlitz- 
Thouless-Berezinskii (KTB) transition in an overdamped Josephson junction array. At temperatures 
slightly above the KTB transition temperature, the noise is white for f < f$ and scales as 1/f for 
/ > /{. Here oc £ -z , where £ is the correlation length and z is the dynamic exponent. Moreover, 
when all frequencies are scaled by /j, data for different temperatures and frequencies collapse on to 
a single curve. 

PACS numbers: 74.50. +r, 74.40. +k 
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Arrays of Josephson junctions have been used exten- 
sively as a model system to study the effects of order pa- 
rameter phase fluctuations on the superconducting tran- 
sition in two dimensions. Such arrays can be fabricated 
with a high degree of uniformity and their relevant pa- 
rameters can be accurately determined. It is widely ac- 
cepted that the zero field transition is described by the 
Kosterlitz-Thouless- Berezinskii (KTB) theory Q-|| and 
its extension to non-zero frequency [0-[7|. According to 
this theory phase coherence is established throughout 
the sample below a temperature Tktb, and the system 
is superconducting. For temperatures above Tktb but 
below the bulk transition temperature, even though indi- 
vidual islands are superconducting the array is not. The 
thermal excitations, vortices and antivortices, that trig- 
ger this phase transition are topological defects of the 
order parameter. Below Tktb vortices and antivortices 
bind in pairs to produce a vortex dielectric, while above 
Tktb the pairs dissociate to form a vortex plasma. In 
the vortex plasma phase, one can identify a character- 
istic length, £, as the average separation between free 
vortices; as T — > Tj^ TB ,^ diverges. Thermal fluctuations 
that perturb the vortex density away from its equilib- 
rium value relax through some local dynamic process. 
Thus, associated with the characteristic length £ there is 
a characteristic time r (or an inverse characteristic fre- 
quency /jT 1 oc r) corresponding to the time required for 
the disturbance to propagate across the distance £. As 
£ diverges, so does r (f^ 1 ), signifying critical slowing 
down. In general r = t {^/^ ) z , where the exponent z 
depends on the dynamics of the relaxation and tq and £o 
are non-universal time and length scales characteristic to 
the specific sample. For simple diffusion, z = 2. 

Previous experimental studies have involved both elec- 
trical resistance and two-coil mutual inductance 
techniques, both of which apply an external force 
to the system and are generally confined to a specific 
frequency. Because the transition to the resistive state is 
determined by the dissociation of vortex pairs by thermal 
fluctuations, these external forces, which also dissociate 
pairs, affect one's ability to study the intrinsic critical 
phnomena near the true thermodynamic transition tem- 
perature. In this paper we employ a non-invasive probe 
to study the transition in equilibrium. Specifically, we 
measure the spectral density of magnetic flux noise [18| , 
S$(f), over a frequency range of more than five decades. 
We find that S$(f) is white for f < and scales as 1/f 
for / > /f . In addition, by plotting fS$(f) versus ///t 
we show that the data collapse in a manner consistent 
with dynamic scaling. 

The 1 mm x I mm array consists of 0.2 /jm-thick 
cross-shaped niobium islands patterned on top of a 0.3 
jtim-thick copper film [Fig. 0(a)]. The islands form a 
square array with a lattice constant of 10 /xm, and the 
junctions are 4 /xm wide and 2 /im long. The mea- 
surement apparatus, originally used to study vortex mo- 



tion in high-temperature superconductors |20j], involves 
a Nb-based Superconducting QUantum Interference De- 
vice (SQUID) attached to a cold stage inside a vacuum 
can surrounded by liquid 4 He. A mu-metal shield reduces 
the static magnetic field to less than I /iT. The array, 
equipped with current and voltage leads, is mounted a 
distance d (< 100/im) away on a variable temperature 
stage. The SQUID is a square washer with inner and 
outer dimensions ii = 180 ^m and i Q = 900 //m [Fig. 
0(b)], and is operated in a flux-locked loop. The out- 
put signal is proportional to the change in magnetic flux 
through the SQUID induced by vortex motion in the ar- 
ray. 

The inset to Fig. [2] shows the differential resistance, 
dV/dl, of the array versus temperature, T, measured at 
zero bias current with an rms current of 10 fj,A at a fre- 
quency of 47 Hz. The initial drop in the resistance at 
approximately T — 8 K marks the bulk transition tem- 
perature of the niobium islands. As the temperature is 
lowered, a resistive plateau develops and is followed by 
a second precipitous drop, which is the KTB transition. 
These data are similar to those obtained in previous ex- 
periments §-0. To obtain the average critical current 
per junction, « C (T), we measure dV/dl as a function of 
the static bias current. We take the current at which 
dV/dl is a maximum as Ni c (T), where N is the number 
of junctions across the width of the array [ pTp^ . 

Figure]^ shows the spectral density of flux noise, S$(/), 
for 15 temperatures above Tktb- At each temperature 
5«(/) is white for / < f 6 {T) and cx 1/f for / > / C (T). 
We define fe(T) as the intersection of lines through 
the white and 1/f noise regions as shown in Fig. |^ for 
T = 1.978 K. Qualitatively, we understand the differ- 
ence between the two frequency regimes as follows. For 
/ > we probe the system at a time scale shorter than 
that required for vortex density disturbances to travel 
over the correlation length. Thus, the system appears to 
be critical. On the other hand, for / < the time scale 




FIG. 1. (a) Schematic of array. Crosses are niobium is- 
lands; area between crosses is copper film, (b) Schematic 
representation of SQUID, with inner and outer dimensions 
ii and i , a distance d from array; dashed square has area 
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FIG. 2. Spectral density of magnetic flux noise, S$(/), ver- 
sus frequency for 15 temperatures above Tktb; scatter at 
higher temperatures is due to subtraction of SQUID noise. 
Dashed lines have slope —1 and 0. Inset shows dV/dl versus 
T. 

is longer than the characteristic relaxation time, so that 
the system appears disordered. In this picture, the dis- 
tinctions between the 1// (critical) and white (disor- 
dered) power spectra reflect the fundamental differences 
in the relaxation dynamics of these two very different 
thermodynamic states. It is worthwhile to point out 
that theories based on the two-dimensional XY model 
and time-dependent-Ginzburg-Landau dynamics predict 
a l// 2 behavior for f > f$ |2^;|2^]. This is in marked 
contrast to the 1// behavior observed in our experiment. 

To interpret our data, we now present a brief discus- 
sion of a scaling theory for the flux noise measurements. 
In our geometry [Fig. |l|(b)], h is the direction perpen- 
dicular to the plane (h = 0) defined by the array. We 
denote a three-dimensional vector by (x, h) , where x is 
the component of the vector in the plane of the array. In 
the absence of the SQUID, the perpendicular component 
of the magnetic field, B±, induced at point (x, h) at time 
t by a vortex density distribution p v (x' , i) at h = is 

Bj_(x,h;t) = $ d 2 x'K(x- x",h)p v (x',t). (1) 



Here $ = hc/2e and K(x-x',h) cx h/(\x- x'\ 2 + h 2 f' 2 
for \x — x'\ 2 + h 2 ^> a 2 [^3), where a is the lattice con- 
stant of the array. The total flux, detected by the 
SQUID is $(£) = J d 2 xB J _(x 1 h = d;t), where the inte- 
gral is performed over the effective area of the SQUID, 
l 2 s = £i£ . In fact, magnetic field lines produced by vor- 
tices near the SQUID are distorted by the presence of the 
superconducting washer, modifying the form of K. How- 
ever, we do not expect this modification to change the 
scaling dimension of K from two (that is, if we replace 



x,x\d by x/C d/£, we still expect K(\x — a? \,d) to 
become £,~ 2 K((x — x')/^, d/Q so that our scaling argu- 
ments remain valid. Thus the flux-flux correlation func- 
tion, C*(t) = ($(t)$(0)>, is 

C 9 {t) = $ J d 2 xd 2 yd 2 x'd 2 y'K(x - x' , d) 

xK(y-g,d)( Pv (x',t)p v (y',0)). (2) 

Here, ( ) denotes the thermodynamic average, and the 
unprimcd integrals are taken over the effective area of 
the SQUID. Next, we assume a scaling ansatz for the 
vortex density-density correlation function: 



( Pv (x\ tK(^, 0)) - C'F^t/T, \x' - L/C). 



(3) 



Here we have used the fact that the scaling dimension 
of the vortex density is two at the KTB transition (that 
is, p v oc £~ 2 ). In Eq. (||), F\ is a scaling function and 
L is the smaller dimension of the array. For Josephson 
junction arrays, 



t = £oexp (ty\/ T '- T KTB 
T' = k B T/E,,(T) = 2ek B T/M c (T), 



(4a) 
(4b) 



and Ej(T) is the Josephson coupling energy per junction 
Q. If we substitute the scaling ansatz, Eq. (||), into 
Eq. (||) and perform proper rescaling of the integration 
variables, one can show that 



C*(t) = $ 2 F 2 (t/T,d/£,L/£,4ff/U 



(5) 



where F% is a new scaling function. The noise spectrum, 
S*(f), is defined as 5$(/) = / dt exp(i27r/i)C$(i). Sub- 
stituting Eq. da) in this expression, we find 



/5*(T) - ^ 2 F(f/f^d/CL/C£es/0- 



(6) 



Here F is another scaling function, the form of which we 
determine from a subsequent data collapse. 

If we ignore the dependence on d/C L/C and l e ff/£, 
Eq. (0) predicts 



fS 9 (T) = *gf ((///„) cxp(6z/ v /T' - T^ TB )) , 



(7) 



where we have explicitly put in the temperature de- 
pendence of /j. In Eq. (j?]), Tktb and bz are un- 
known. To determine T", we fit i c (T) between T = 
1.52 K and T = 1.93 K using the expected tempera- 
ture dependence i c (T) = i c (0) exp(— aVT) ^5|, find- 
ing i c (0) = 0.133 A and a = 9.14 KT 1 / 2 , and extrap- 
olate to higher temperatures. Equation (|t]) predicts that 
by choosing the correct values for Tktb and bz and 
plotting fS$(f) versus the scaled frequency f/rofe = 
fexpfbz/y/T' — T^ TB ), we should obtain a collapse of 



3 



the raw data onto a single curve. Our procedure for 
data collapsing is to make iterative changes in bz and 
Tktb until the best data collapse is obtained. The final 
result is shown in Fig. where the fitting parameters 
were bz — 4 and Tktb = 1-63 K. The quality of the 
data collapse is strongly affected by the choice of Tktb 
but relatively weakly by the choice of bz. One can com- 
pensate an increase (decrease) of Tktb of a few mK by a 
decrease (increase) in bz of about 10%. However, changes 
in Tktb beyond a few mK or changes in bz of more than 
about 10% always result in a lower quality data collapse. 
We can obtain a more accurate estimate of bz through 
the temperature dependence of /f(T). Using Eqs. (Q) 
and fe(T) — /o[£(T)/£ ]~ 2 , we expect a straight line of 
slope —bz on a plot of ln(/^) versus 1/yT — T KTB . In- 
set (a) to Fig. ||, using Tktb = 1-63 K, shows the re- 
sultant straight line which yields bz — 4.10 ± 0.04 and 
fo = (2.1 ±0.3) x 10 6 Hz. 

We note that although /o is in the megahertz frequency 
range, the exponential dependence of on T' — T KTB im- 
plies that is reduced to the order of 1 Hz even when 
T is about 0.2 K above Tktb- This sensitivity is illus- 
trated by the fact that the apparent transition temper- 
ature, Tktb (47 Hz), associated with the sharp rise in 
dV/dl measured at 47 Hz (inset, Fig. ||) is substantially 
higher than the static value Tktb = 1-63 K inferred from 
our noise data. However, these values are quite con- 
sistent since Tktb (47 Hz) is the temperature at which 

= 47 Hz. Using the above values of Tktb, /oj an d bz, 
we predict Tktb (47 Hz) = 1.94 K, in good agreement 
with the dV/dl data (inset, Fig. ||). 

In addition to finding bz , we can extract the critical 
exponent z in the following manner. At low frequencies, 
a two-dimensional sample with a linear conductivity a 
produces white noise with spectral density oc ksTa. 
In the KTFJ regime for / <C /f , o is proportional to 
(Pv) — £ 2 > reflecting the intuitive idea that the low 
frequency conductivity is inversely proportional to the 
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FIG. 3. fS<f>(f) versus f/rofe for flux noise shown in 
Fig. § Tktb = 1.63 K, bz = 4. Insets show ln(/ e ) versus 
l/^/T' - T^ TB with line of slope -4.10 and ln(S* /T) versus 
l/<s/T' - T£ TB with line of slope 4.13. 



density of free vortices. Therefore, for / <C 

S3? « k B Te exp(26/^/T'-T^ TB ) . (8) 

To test Eq. (||), we define S% as the horizontal line 
drawn through the low frequency data at each temper- 
ature, as shown for T = 1.978 K in Fig. E| and plot 
ln(S|/T) versus 1/VT'-T KTB . Inset (b) to Fig. |, 
using Tktb = 1-63 K, shows the resultant straight line 
which yields 2b = 4.13±0.04. Combining this result with 
the value bz = 4.10 ± 0.04 from the temperature depen- 
dence of /{ yields the dynamic exponent z = 1.98 ±0.03. 

In summary, we have used measurements of magnetic 
flux noise to study the equilibrium KTB transition in an 
array of overdamped Josephson junctions. We emphasize 
that confirmation of the critical dynamics of the KTB 
transition requires measurement of frequency and tem- 
perature dependent properties such as those presented 
above. Using bz — 4 and Tktb = 1-63 K, we have shown 
that the data collapse over more than five decades in 
frequency, confirming the predictions of dynamic scal- 
ing. In addition, from the temperature dependences of 
/j and S% , we have extracted the dynamic exponent z — 
1.98±0.03. Our experimental finding S$(/) cx 1// in the 
critical frequency regime is inconsistent with theoretical 
predictions based on time-dependent Ginzburg-Landau 
(TDGL) dynamics with the classical two-dimensional XY 
model p3| , p4| ]. The failure of this model is surprising 
since TDGL dynamics does predict z = 2, in agree- 
ment with our extracted value. It is unclear to us why 
£*(/, T) is independent of d/£, L/£, and l e fi/£- A possi- 
ble cause may be that the data in Fig. are in the regime 
£(T) < £ e ff < L. Our evidence for this restriction is that 
at temperatures between Tktb and 1.825 K (the lowest 
temperature for which data are shown in Fig. |J), we ob- 
serve discrete jumps in the flux threading the SQUID. 
We interpret this behavior as the motion of a single vor- 
tex under the SQUID, implying that £(T) > £ e g, and 
conclude that £(T) < £ g for the temperatures referred 
to in Fig. |. If we assume £(T = 1.825 K) < £ eS , we 
deduce fo < 0.2 ^um [Eqs. (Q)], a value that is consider- 
ably smaller than the commonly accepted lore fo = a - 
Setting fo < 0.2 /xm, we deduce f < 2 /im for the high- 
est temperature (2.379 K) data shown in Fig. ||. The 
fact that scaling persists to this temperaure where the 
correlation length is smaller than the lattice constant a 
is somewhat disturbing. Further work is needed to re- 
solve this issue, as well as the lack of dependence of F 
on d/£, T/f , and i e s/^, and the 1// behavior of the flux 
noise in the critical regime. 

In closing, we note that we have carried out similar 
measurements on two other arrays, one of square ge- 
ometry and the other of triangular geometry, using the 
SQUID described above as well as two other SQUIDS of 
different geometry. When the apparatus was cooled in a 
magnetic field (< 1 jtiT), the data were similar to those 
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described above with minor differences in the high fre- 
quency behavior. When a small field was applied, the 
noise data exhibited very different behavior. These ad- 
ditional measurements will be described elsewhere. 
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